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Abstract 

A  new  numerical  simulation  method  using  the  finite  element  methodology  (FEM)  is  presented  to 
study  electromagnetic  scattering  due  to  an  arbitrarily  shaped  material  body  doped  randomely  with  thin  and 
short  metallic  wires.  The  FEM  approach  described  in  many  standard  text  books  [1,2]  is  appropriately 
modified  to  account  for  the  presence  of  thin  and  short  metallic  wires  distributed  randomly  inside  an  arbitri- 
ally  shaped  material  body.  Using  this  modified  FEM  approach,  the  electromagnetic  scattering  due  to 
cylindrical,  spherical  material  body  doped  randomly  with  thin  metallic  waires  is  studied. 

I.  Introduction 

Electromagnetic  (EM)  scattering  properties  such  as  monostatic/bistatic  radar  cross  section  of  a  homo¬ 
geneous  material  object  can  be  controlled  or  modified  by  embedding  metallic/non-metallic  inclusions  in 
the  object.  These  inclusions  can  be  in  the  form  of  short  metallic  wires,  small  thin  metallic  plates,  or  small 
metallic  particles  of  various  shapes  [3].  The  EM  scattering  properties  can  be  controlled  by  selecting  proper 
size  and  shape  of  these  inclusions.  Recent  advances  in  fabriction  technology  have  allowed  for  the  inclu¬ 
sion  of  nano-scale  metallic  wires  in  these  objects.  However,  current  fabrication  techniques  used  have  no 
control  over  the  spacing  and  arrangement  of  these  small  sized  wires  in  a  binding  medium.  As  a  result, 
these  short  wires  are  arranged  in  a  arbitrary  fashion.  It  is  important  and  economically  advantageous  to 
know  the  EM  scattering  properties  of  these  objects  populated  randomly  with  these  wires/particles  prior  to 
its  actual  fabrication.  In  this  work  we  describe  a  modifed  FEM  approach  to  determine  the  EM  scattering 
properties  of  randomly  dispersed  thin  and  short  metallic  wires  in  an  arbitraily  shaped  material  object. 

II.  Numerical  Modeling 

In  this  section,  first,  we  present  a  brief  outline  of  the  FEM  approach  to  estimate  the  EM  scattering  from  an 
arbitrailly  shaped  material  object  without  any  metallic  wires  or  other  inclusions.  Then,  we  explain  the 
steps  to  modify  this  simple  FEM  approach  to  account  for  presence  of  randomly  populated  thin  and  short 
metallic  wires  in  the  object.  Figure  1  shows  an  arbitrily  shaped  material  object  of  permittivity  £r  and  per¬ 
meability  |Ir  populated  randomly  with  small/nano  size  metallic  wires.  To  estimate  the  EM  scattering  from 
the  object  shown  in  Figure  1  without  any  metallic  wires,  the  entire  scattering  volume  is  split  into  two 
regions  by  enclosing  the  object  by  a  ficticous  surface  as  shown  in  Figure  1.  The  electric  field  inside 

the  surface  ,S'|  satisfy  the  vector  wave  equation 

Vx(l/| \ir)VxE  -  k20erE  =  0  (1) 

Following  the  usual  steps  invovled  in  the  FEM  formulation,  the  equation  (1)  can  be  reduced  to  [4] 

J (( 1/|H JVxT'VxE1 -k20£rE  *T)dv  =  -$T*  nx((l/yLr)VxE)dsl  (2) 

V  5] 

where  T  is  a  vector  testing  function  and  n  is  the  unit  outward  normal  to  the  surface  ,S'|  .  If  you  introduce 

the  surface  current  J  as  an  additional  unknown,  then  the  right  hand  side  of  equation  (2)  can  be  simplified 
to 
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Figure  1  Geometry  of  arbitrary  shaped  material  object  embedded  randomly  with 
thin  and  short  metallic  wires. 

The  equation  (2)  can  then  be  written  as 

J((l/[ir)Vxf*  VxE  -k20£rE  •T)dv-j(n\i0$T*Jdsl  =  0 


(3) 


V 


S 1 


In  equation  (3),  in  addition  to  the  unknown  electric  field,  the  surface  current  J  is  introduced  as  an  extra 
unknown  and  to  generate  additional  equations,  continuity  of  tangential  electric  field  across  the  fictious  sur¬ 
face  boundary  can  be  used.  Hence, 


E 


■)  -) 

=  -  VxF-j'co|i0A  +  l/(/(0£0)VV»A  +  Einc 


(4) 


where  A  and  F  are  the  vector  potential  functions  for  the  region  outside  the  surface  Sl  and  Einc  is 
the  incident  field.  To  facilitate  the  solution  of  equations  (3)  and  (4)  the  region  inside  the  ficticious 
surfae  3,  is  discretized  into  terahedron  and  the  electric  field  over  each  tetrahedron  can  be 


expresse 


N 


-> n 


E(x,  y,  z)  =  Yj  H  bmWm(x’y>z) 


n  =  1  m  =  1 


(5) 


n  n 

where  b  and  Wm(x,  y,  z.)  arc,  respectively,  the  unknown  amplitude  of  electric  field  and  the  vector  basis 

th  tli  $  -A 

function  associated  with  the  m  edge  of  n  tetrahedron.  The  surface  currents  J  and  M  required  for  the 

calculations  of  vector  potentials  can  also  be  expressed  in  terms  of  vector  basis  functions  W  as 

3  3 

2  i,dt  x  » 
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j=  YIinx  Wi> 
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(6) 


i  =  1 
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Using  the  expressions  (5),  (6),  and  the  Method  of  Moments  (MoM),  the  equations  (3)  and  (4)  arc  con¬ 
verted  into  set  of  simultaneous  equations: 

S 1 1_  _^12_  \b\  _  [0 

^21  $22_  _[{]_  _[V] 

If  N ed  are  the  number  of  tetrahedron  edges  and  N  t  are  the  edges  lying  on  the  surface  ,  then  the  coef- 

n 

ficient  matrix  is  of  the  order  (Ned  +  Napt)  X  (Ned  +  Napt)  .which  can  be  solved  for  bm  and  /■ .  From  the 

knowledge  of  the  surface  currents,  the  EM  scattering  properties  of  an  arbitrarilly  shaped  object  are  esti¬ 
mated. 

The  simulation  technique  presented  so  far  estimates  the  EM  scattering  from  a  homogenous  material 
object.  To  estimate  the  EM  scattering  due  to  an  arbitrarilly  shaped  body  populated  randomly  with  short 
thin  wires  we  follow  the  procedure  decribed  below.  Let  us  assume  that  when  the  object  is  embedded  with 
thin  wires,  a  wire  occupies  one  of  the  positions  of  edges  out  of  total  N  ed  edges.  To  populate  randomley 

the  homogeneous  object  with  Nwjre  number  of  wires,  we  generate  Nwire  random  numbers  between 
{ 1  ~  N ed } .  The  edges  corresponding  to  these  random  numbers  are  assumed  to  be  occupied  by  thin 

n 

metallic  wires.  Consequently,  bm  coefficients  on  these  edges  will  be  zero.  The  modified  matrix  equation 

for  the  object  randomly  populated  with  thin  metallic  wires  can  be  simply  obtained  by  eliminating  rows 
and  columns  corresponding  to  those  edges  where  metallic  wires  arc  assumed  to  be  present.  Hence  the 
matrix  equation  (6)  when  applied  to  the  object  loaded  with  thin  metallic  wires  gets  modified  to 

Srecl_  [p\  =  vred_  ^ 

where  the  Sred  and  vred  arc  the  matrices  obtained  from  and  |~JJ  by  eliminating  the  rows  and 

columns  corresponding  to  the  edges  where  metallic  wires  arc  assumed  to  be  present. 

III.  Numerical  Results  and  Discussion 

For  a  numerical  experiment  we  consider  a  material  sphere  of  radius  k{)a  =  1.0  with  permittivity 
Er  =  4.0  - j'0.0  ,  permeability  |l  =  l.O-j’O.O  illuminated  by  a  plane  EM  wave.  To  estimate  bistatic 
radar  cross  section  of  the  material  sphere  we  assume  that  the  sphere  is  illuminated  by  a  plane  wave  with 

incident  angle  0;;!  =  180°,  (j)  M  =  0°.  To  facilitate  the  bistatic  RCS  calculation  the  sphere  is  decretized 

using  the  COSMOS/Geostar  as  is  shown  in  Figure  2(a).  Using  the  procedure  described  above  the  bistatic 
RCS  of  the  material  sphere  is  calculated  and  presented  in  Figure  2(b).  The  number  of  tetrahedron  used  to 
discretized  the  sphere  were  123  resulting  in  Ned  =  238  and  N  t  =  132 .  The  numerical  data  presented 

in  Figure  2(b)  arc  confirmed  with  the  other  published  results. The  agreement  between  the  results  obtained 
using  the  present  procedure  and  other  published  results  confirms  the  validity  of  the  present  formulation. 
Now  to  estimate  the  bistatic  RCS  for  the  material  sphere  doped  with  Nwjre  number  of  very  thin  metallic 

wires,  Nwire  random  numbers  uniformly  distributed  between  1  and  N ed  =  238  arc  generated.  If  P d 
is  the  percentage  of  doping,  then  the  number  wires  to  be  used  for  doping  can  be  calculated  from 


Figure  3(a)  Material  sphere  (  k^a  —  1 ) 
populated  randomly  with  thin  metallic  wires 
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Figure  3(b)  Bistatic  RCS  of  material  sphere  shown 
in  Figure  3(a) 
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For  the  present  example,  with 


P d  =  5  and  N ed  =  238  choices  the  number  of 


wires  used  for  doping  were  Nwjre  =  12  .  The  edges,  corresponding  to  these  12  random  numbers  (uni¬ 
formly  distributed  between  1  and  238  )  are  assumed  to  be  thin  metallic  wires.  The  doped  material  sphere 
with  the  randomely  oriented  wires  is  shown  in  Figure  3(a).  Following  the  FEM  procedure  described 
above,  the  bistatic  RCS  for  the  new  structure  is  estimated  and  is  shown  in  Figures  3(b).  Various  numer¬ 
ical  examples  were  simulated  for  varying  concentrations  of  wires,  their  orientations  and  placing.  These 
results  will  be  discussed  at  the  time  of  presentation. 

IV  Conclusions 

A  new  and  modified  finite  element  methodology  (FEM)  has  been  successfully  developed  to  study 
electromagnetic  scattering  from  a  arbitrarily  shaped  material  object  populated  randomly  with  thin  metallic 
wires.  Bistatic  RCS  of  material  objects  in  spherical,  cylindrical,  and  cubical  shapes  randomly  populated 
with  thin  metallic  wires  have  been  studied. 
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